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Abstract 


Many phenomena in various fields of physics are simulated by parabolic 
partial differential equations with the nonlocal initial conditions, while there 
are few numerical methods for solving these problems. In this paper, the 
Ritz—Galerkin method with a new approach is proposed to give the exact and 
approximate product solution of a parabolic equation with the nonstandard 
initial conditions. For this purpose, at first, we introduce a function called 
satisfier function, which satisfies all the initial and boundary conditions. The 
uniqueness of the satisfier function and its relation to the exact solution are 
discussed. Then the Ritz—Galerkin method with satisfier function is used to 
simplify the parabolic partial differential equations to the solution of algebraic 
equations. Error analysis is worked by using the property of interpolation. 
The comparisons of the obtained results with the results of other methods 
show more accuracy in the presented technique. 
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1 Introduction 


Various problems arising in geology [19], heat conduction [6, 7, 20], plasma 
physics [16], chemical engineering [12], thermostatic [30], and hydrodynam- 
ics [11], can be reduced to nonclassical initial-boundary value problems. The 
investigation of a nonclassical problem with nonlocal initial condition is con- 


sidered in this article. Nonclassical problems with nonlocal initial conditions, 
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which are generalizations of the classical or time-periodic problems, can be 
used in science with better results than the classical initial condition. Nonlo- 
cal initial conditions are useful in the modeling of radionuclides propagation 
in Stokes fluid, sewage causing pollution processes in rivers and seas, and 
diffusion in porous media; see [17, 25, 29]. 

The following parabolic equation with nonlocal initial condition is con- 
sidered in this paper: 


Ov _ 02u 
Ot Oy? 


with boundary conditions 


+¢(7,t), O<y<1, O0<t<T, (1) 


v(0,t) =Ao(t), O<t<T, (2) 
v(1,t)=h(t), O<t<T, (3) 
and the nonlocal initial condition 
N 
u(y,0) = S° Aj (yu(7, 73) + x(); 
j=l 
O<y<1l, 0<7, <TIhn<:::-<Ty=T, (4) 


The function v should be found, but ¢, ho, hi, x, and A; are known functions. 

The multiplier method, semi-group theory, the maximum principle, and 
potential theoretical representation of the solution are applied to prove the 
existence and uniqueness of the above problem; see [8, 10, 3, 4, 5, 9, 24]. 
There are many publications about numerical methods for parabolic par- 
tial differential equations with classical initial condition and with nonlocal 
boundary conditions, but there are few numerical methods for nonstandard 
initial condition problem presented here [13, 14, 15, 22]. Therefore, it is sig- 
nificant to propose numerical methods for this problem. Dehghan [13, 14] 
concentrated on finite difference schemes to solve problem (1)—(4). The im- 
plicit collocation technique is applied for solving the above problem in [15]. 
Recently, in [22], the authors analyzed finite difference schemes for solving 
(1)—(4). 

Several partial differential equations are numerically solved by Ritz and 
Galerkin method [31, 32]. In the Ritz-Galerkin method, applying the ap- 
propriate satisfier function has been important; see, for instance, [33, 34, 35, 
2, 21, 26, 1]. Utilizing satisfier function that fulfills all the problem condi- 
tions helps to satisfy the problem conditions quickly, in addition, it leads to 
a smaller system of algebraic equations and reduces the time of computation. 

In this paper, we concentrate on satisfier function of (2)-(4). The unique- 
ness of the satisfier function and its relation to the exact solution are discussed 
in some theorems. Moreover, the Ritz-Galerkin method in the Bernstein 
polynomial basis applying the satisfier function is utilized to give an ap- 
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proximate solution of (1)-(4). The comparison of the obtained results with 
those results obtained by [13, 14, 22] shows more efficiency of the presented 
technique. 

This paper is separated into the following sections: In Section 2, the 
satisfier functions for (2)—(4) are introduced. Theorems on the uniqueness of 
product satisfier functions are proven and used to obtain the exact solution 
of the problem. Moreover, the Ritz—Galerkin method with satisfier function 
is applied for solving the problem numerically. Some results concerning the 
error analysis are obtained in Section 3. In Section 4, we apply these results 
to solve three nonlocal initial condition problems. Section 5 is dedicated to 
the conclusion. 


2 Satisfier function for Ritz—Galerkin method 


The solution v(7,t) is approximated with the following series in the Ritz 
method: 


n m 
B(7,t) = S> Yd bizoug (7, t) + n(7, t), (7, #) € [0,1] x (0,7), (5) 
i=0 j=0 
in which o;;(7,t) = y¥v - 1) I(t — T;)pi(y)e;(t), where p;(y)p7(t) are 
basis functions. The function nty,t) is called the satisfier function. In this 
study, we take p;(7) the Bernstein polynomials in [0, 1]. What is important in 
the Ritz method is to find a proper satisfier function. The satisfier function 
is an arbitrary function that satisfies all problem conditions [35]. Satisfier 
functions are not unique, and the construction of these functions is often 
based on interpolation techniques. Here we focus on the construction of 
satisfier function for one-dimensional parabolic equation with the nonlocal 
initial condition. 
We assume that the compatibility conditions 


N 

ho(0) = x(0) + 5° A; (0)ho(T5), (6) 
j=l 
N 

hi (0) = x(1) + Dy dj (1) hi (Tj), (7) 


are satisfied and define the satisfier function 7(7, t) that satisfies the boundary 
conditions (2), (3), and the nonlocal time weighting initial condition (4) when 


N 
ho(0) — S>Ai(qho(Zj) #0, OS 7<1, 
j=l 
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as 
(7, t) = a(y)ho(t) + yhi(t), 
when 
N 
)- So A(Ym(T;) 40, OS <1, 
j=l 
as 
ny, t) = (1 — y)ho(t) + b(y)hi (4) 
and when 


N N 
ho(0) — }>As(a)ho(T3) = ha(0) — AAT) =0, OS 7S 1, 


n(y,t) = (1 — y)ho(t) + rit) + c(t)x(y); 


_ x7) ~ 20 ) + pe As (ha (Ty) 
ho(0) — 37h Ay (Who (Tj) 


x(y) — (1 = 7) ho (0) + 35hy Ay (7)(1 — ho (Ty) 


ans nO Ep yma (T) 
N 


Theorem 1. If at least one of the values ho(0) — eu Aj (y)ho(Z;) and 
hy(0) — pa 1A;(7)hi(T;) are not equal to zero in [0,1], then the product 
satisfier funetion for (2) — (4) is unique. 

Proof. Suppose that 


N 


hy (0) — $5 As(y)ha (Ty) # 0 (8) 


j=l 


and that the product functions A;(y)Ae(t) and Fi (7) F2(t) satisfy (2) — (4). 
Then 


A, (1)Ao(t) = Fi(1) a(t) (9) 
Aj (7) A2(0) — Fi(y)F2(0) = DLAs (Ar (7) A273) — Fi(y)Fo(7;)). (10) 


From (9), in particular, we have 
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Ay (1)A2(Tj) = Fi(1)Fa(2j), Ar (1)A2(0) = Fi (1)£2(0). (11) 
By multiplying (9) by (10) and using (9) and (11), we have 


(Ai (7) Aa(t)—Fi (7) Fa(t))ha (0) = (Ar (7) Aa(t)—Fi (7) Fa) (95 Ag (ha (T5))- 
Thus 


N 
(Ai (1) Ao(t) — Fi (7) Fa(t))(h (0) — (95 g(a) ha (Ty) = 0; 


therefore from (8), we get 


Ai (y)Aa(t) = Fi(y) Fat). 


In a similar way, the other case ho(0) — . Aj (y)ho(T;) 4 0 is proved. 


Theorem 2. Suppose that hi(0) — ae Aj(y)hi(Z;) #0, O<y<1 and 
that 


Me 


N 
[hi (0) — )_ Aj (0)hi(T;)]Ro(t) = [ho(0) — y Aj (O)ho(T;)]hi(t). (12) 


j=l 


Then the unique separable satisfier function for (2) — (4) is given by 


x(a)ha(t) | 
hy (0) — 77h Ag (aD) 


n(7; t) = 


Proof. From compatibility condition (7), the function (7,t) satisfies the 
boundary condition (3) obviously. Now let ho(0) — wy Aj (O)ho(T)) F 0. 
Then by using (12) and the compatibility condition (6), 7(y,t) satisfies the 
boundary condition (2) and when ho(0) — pe A; (0)ho(T;) = 0, from (12), 
we have ho(t) = 0. Thus from the compatibility condition (6), we obtain 
x(0) = 0 and condition (2) also is fulfilled. In addition, n(7,t) satisfies the 


initial condition (4), because 
_ xin) 
ha (0) — Dy AV (T}) 
_ x(r)hx(0) — x7) Sen As (Wha (Ty) + x1) pea Ag (Yaa (Tj) 
hy (0) — phy Ay (ha (Ty) . 


n(7; 0) 


Hence 
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Theorem 3. Suppose that ho(0) — ee Aj (y)ho(Z;) 4 0 and that 


N 


N 
[ha (0) — $7 Aj (Ox (T;)]ho(t) = [ho (0) — S A4(0)ho (Ty)] 1 (4). 


j=l j=l 
Then the unique separable satisfier function for (2) — (3) and (4) is given by 


: x(7)ho(t) 
my t) = ho(0) — Djs As(*)ho(Ti) 


Proof. The proof is similar to the previous theorem. 


Now we consider the approximation (5) and apply the following Galerkin 
equations to obtaining coefficients b;;: 


0b 076 ; : 
< BE Dye PVD» Pisa MP5. n(O) P= 0, (i - 0,...,7), (j — 0,...,m), 
(13) 
where <,> is defined as 
0b 86 
& at _ a7? _ #7, t), Pin (Y)P3,m(t) > 


1? ao 8 
= ff He — oO Mri psom te 


The coefficients b;;(¢ = 0,...,n)(j = 0,...,m) are obtained with solving 
algebraic equations system (13). 

Now let ho(t) = hi(t) = 0 and let a product form for the satisfier function 
be taken. Then it can be denoted as 


n(y,t) = f(y)g(t). (14) 


Substituting the function (14) in to conditions (2)-(3), gives the following 
conditions on f(7): 


f(0)=0, f(1)=0. 
On the other hand, from the compatibility conditions (6) and (7), we have 
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Therefore, we can conclude that one of the options for the function f(y) can 
be x(y), and for fulfilling the initial condition (4) we choose g(t) such that 


N 
9(0) =1+ D/As()9(T)- (15) 
Therefore 
n(7,t) = x(y)9(t), (16) 


with condition (15) is a satisfier function for (2)—(4). 


It is noteworthy that if ho(t) = hi(t) = 0 and x(y) #4 0, then we can 
approximate the product solution as 


(7, t) = x(y)h(t) = x(7) > Pipin(t). (17) 


This approximation supplies higher adaptability in the nonlinear initial and 
boundary conditions. Then, with the Galerkin equations and 


n N n 
> b; Pin (0) ~~ (1 + > Aj (y) » b:0i,n(T;)) = 0, (18) 
i=0 j=l i=0 


the coefficients b; are specified. If we choose Bernstein polynomials for basis 
functions p;n(t), then from properties of Bernstein polynomials, equation 
(18) is equal to the following equation: 


N 
bo — (1+ bn S- Aj(7)) = 0. 


On the other hand, when X(7} is a constant, the separable satisfier function 


is also the solution of the problem(1)—(4) and the exact solution will also be 
obtained from (16) with 


x 


t o(7, 8) exp(—X} 5) 
g(t) = exp Saf ex as +0); 
x(7) “Jo x(7) 
such that c can be obtained from the following condition: 
N 
9(0) =1+ >) As(y)9(T}). (19) 
j=l 


Also if ho(t) = hi(t) = x(y) = 0, then we can try to find the following 
approximation for the solution of the problem: 
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O(y, t) = A(y)g(t) = Aly) 95 dipin(®), 
i=0 


and A(7) is selected such that A(0) = A(1) = 0. For example, y(y — 1) and 
sin(7y) are two selections for A(y). Note that some times ¢(7,t) in (1) can 
help us to obtain an appropriate selection for A(y). Then, the expansion 
coefficients b; are specified with the Galerkin equations and the following 


equation: 
N 


by — bn > Ag(y) = 0. 


j=1 


3 Error analysis 


In the current section, we analyze the error of the proposed method. We 
obtain an estimate of the error norm of the best approximation of a smooth 
function of two variables, on [0,1] x [0,1], by Bernstein polynomials. 


Lemma 1. Let f(y,t) be a sufficiently smooth function in [a,b] x [c,d] 
and let pr—1m-1(y,t) be the interpolating polynomial to f(z,t) at points 
(y:,t;), where y;,, i =0,1,...,n—1, are the roots of the n — 1-degree shifted 
Chebyshev polynomial in [a,b] and t;, 7 = 0,1,...,m—1 are the roots of 
the m — 1-degree shifted Chebyshev polynomial in [c, d]. Then [18, 23] 


(b— a)” a” f(y, t) 
t =, n— m— yt < a 
| (7, t) — Pn—1jm—1(% #) |S oa a ae aan | 
(a— om an Fst) 
a max | ———— | 
mi22m—1 (7,t) € [a,b] x [c,d] ot™ 
(b _ a)"(d _ cy? | OP EOE E (oy: t) | 
min!22(n+m)—2 (+ 4)€ [a,b] x [c,d] Ot™ Oy” : 


Theorem 4. Suppose that fnm(7,t) = yo Dojo bij 045 (7, t) = CTW(y, t) 
is the Bernstein expansion of the real sufficiently smooth function f(y,t) € 
(0, 1] x [0,1]. Then, there exist real numbers c),c2, and c3 such that 


C1 c2 C3 


If —fnm(v t)|l2 < (n+ 1)!22n+1 (m+ 1)!22™4+1 (m +.1)\(n + 1)!22(r+m) +2? (20) 
where ae 
me ~.¢ 
max | cana 82) I< ci, 
(y,t)€ [0,1] x [0,1] Oynrt 
omtl t 
ins | er Lyf) I< ca, 
(y,t)€[0,1] x [0,1] ormri 


oe. fo ee 1 
(y,t)€[0.1}x [0,1] | OLM+TQyntt I= 
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Proof. Let pnim(y,t) be the interpolating polynomial to f at the roots of the 
Chebyshev polynomial. Then from the definition of the best approximation, 
we have 


Ft) — fram (Vs t)\l2 SNF) — Prm (7 tlle: 


Therefore by Lemma 1, we obtain 


| | (F(t) — CP U(y, 1) Payal 


< / / ate eDrnd 


sa io 1 antl f(y, t) 
2) i ee ee 
Jo Jo (m+ 1)!22"+1 (4,4) €[0,1] x [0,1] Oyrtl 
1 an f(y, t) 


(m a 1 )igeeert theo ie [0,1] | otmtl | 
1 ontm+2 fy t) 
tae I|dyat 
(m+ L)!(n + 1)!22+m)+2 (aero ie [0,1] | den +194n41 IP dy 


<| C1 C2 C3 ? 
! (n+ 1)!22P42 7 (mm + 1)122™+1 "(m+ 1)'(n + 1)!22(r+m) +2 


and this completes the proof. 


Remark 1. In the special case when n = m, we have 


C3 ) 1 
(m+ 1122417 (n+ 1)!22-41 


IF Cy, t) — fnn(y t)ll2 < (er + e2 4 


Therefore when considering the upper bound of the error, given by (20), the 
term containing c3 can be neglected and we have 


F(t) — fanl% tlle = Of ). 


1 
(n+ 1)1QnH 


Remark 2. Let f(t) € C"*'[0,1] and let fn(t) = > bipin(t) be the Bern- 
i=0 


stein expansion of f(t). Then 


IFC) — fr(lle = OC ). 


1 
(n aie 1)!22n4+1 
Now let Yn m be the space of bivariate polynomials of degree less than or 
equal to n on y and degree less than or equal to m on t such that the 
Ritz—Galerkin approximation belongs to it. The error norm of the numerical 
results obtained from the Ritz—Galerkin method tends to zero with the same 
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convergence rate as the error norm of the best approximation of the exact 
solution in Y;,,m (27, 28], which is confirmed by the obtained results in the 
next section. Hence we present the following theorem. 


Theorem 5. Suppose that u(y,t) is the exact solution of (1)—(4) and that 
0(y,t) is the approximate solution (5) with (m =n). Then 


leet) 00.1 = OC ayaa) 


4 Numerical examples 


In this section, we report some results of our numerical calculations using the 
methods presented in the previous section for finding the exact and approx- 
imation solutions of (1) — (4). We select Examples 1 and 2 of [13, 14, 15, 22] 
and compare our approximation results. We use the package of Mathematica 
version 11, with the following hardware configuration: Intel(R) Core(TM) 
73 — 2100 CPU, 8 GB of RAM, 64-bit Operating System (Windows 7) for 
numerical calculations in all examples. The approximate L? norm of abso- 
lute error and the maximum norm of the absolute error are, respectively, 
calculated by 


1 oT 1 T 
le, 8) Io = | | e?(,t)dtdy = / / (v(t) — 0(y,2))Pdtdy 


and 
Boe t)| = t) — 6(y,t)]. 
00 jer OO )| eg )=0(,1)| 


Moreover we introduce the notation A, = CHS acee According to the pre- 
vious section, A, has the same asymptotic behavior as the error norm of the 
numerical results, which is confirmed by the results in the tables. 


Example 1. Consider (1)—(4) with N = 2, and let 
o(7, t) = (x? ~ 1) sin(77) exp(—t), 


v(7,0) = vy, T1) — u(y, Ta) + x0), 
xO) =saGyl-e =e), 
ho(t)=0, h(t) = 0, 

T, =0.5, To =1.0, 

Ai(y) = 1, A2(y) = 1. 
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for which the exact solution is [13, 14, 15, 22): 


u(y, ) = sin(ay) exp(—t). 


The results of errors evaluated for different values of n by using the 
method presented in this article with the truncated series (17), together with 
CPU time are shown in Table 1. In Figure 1, the absolute difference between 
the exact and approximate solution for n = 5 is shown. 

In [22], Martin-Vaquero and Sajavicius employed the two-level finite dif- 
ference schemes to solve this problem for different values of h and 7, where h 
and 7 are space and time step sizes, respectively. The best absolute error of 
the approximate solutions at some different points reported in [22] is shown 
in Table 2 and compared with presented method with n = 7. 

The best errors in L? norm obtained in [13] and [14] for solving this 
problem using finite difference schemes with several choices of h, are shown 
in Table 3. 


Table 1: Absolute errors and L? norm of errors for Example 1 


y=t n=1 n=3 n=5 n=7 
0 0 0 0 0 
0.1 8.77648 x 10-2 | 5.50078 x 10-° | 8.63336 x 10-8 | 1.22592 x 10-9 
0.2 1.54731 x 10-! | 2.9701 x 10—-® | 2.10798 x 10-7 | 3.87638 x 10-" 
0.3 1.89538 x 10-1 | 1.10471 x 10-7 | 9.06801 x 10-8 | 5.27107 x 10-% 
0.4 1.88218 x 10-! | 2.3564 x 10-7 | 4.16751 x 10-7 | 1.47424 x 10° 
0.5 1.54818 x 10-! | 2.21861 x 10-7 | 5.25682 x 10-7 | 1.29762 x 10-° 
0.6 1.0049 x 10-! | 6.15521 x 10-° | 8.80609 x 10-9 | 9.66877 x 10- 
0.7 4.12702 x 10-2 | 1.27031 x 10-7 | 4.84358 x 10-7 | 1.08301 x 10-9 
0.8 5.05864 x 10-* | 1.89116 x 10~7 | 2.55435 x 10-7 | 6.65846 x 10-” 
0.9 2.24724 x 10-2 | 7.97571 x 10~® | 1.52571 x 10-7 | 2.07613 x 10-1 
1 0 0 0 0 
Xn 6.25 x 1072 3.25521 x 10-4 | 6.78168 x 107 | 7.56884 x 1071° 
errors in L? norm | 1.3015 x 10-! | 1.35747 x 10-7 | 3.08435 x 10-7 | 9.2949 x 10-9 
CPU Time 0.374 s 0.39 s 0.39 s 0.405 s 


Table 2: Absolute errors for (t = 1) and different values of y for Example 1 


Y Scheme 1 in [22] Scheme 2 in [22] Scheme 3 in [22] Presented method 
0.25) 2.89742 x 10-6 2.32464 x 10-° 9.99528 x 10-14 5.12468 x 10-14 
0.5 | 4.09758 x 10-6 3.28755 x 10-6 1.41368 x 10~1° 7.24995 x 10-14 
0.75) 2.89742 x 10-6 2.32464 x 10-° 9.99524 x 10-14 5.12648 x 10-14 


It is worth pointing out that since the exact solution for u(y,t) is sepa- 


rable, y(y) 4 0, and 


XY 
= 


exp(—t) 


g(t) = 


(l-e-fi —e-hh)’ 


xi is a constant, from (2) and (19), we have 
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_ Error 
| ee x1077 
j = a, 6 


5 


4 


Figure 1: The graph of the absolute error with n = 5 for Example 1. 


Table 3: The best errors in L? norm obtained from different methods and presented 
method for Example 1 


Presented FTCP(13] BTCS[13] Crank-Nicolson[1B|Crandall’s[13] 
method with 
n=6 
8.9 x 109 5.0 x 10° 6.0 x 10° 61x 10° 1.1 x 10-7 
Saulyev I{13] | Saulyev Three-point FTCS | Dufort—Frankel | (1,3, 1) 

II[13] [14] [14] Scheme[14] 
2.0 x 10-7 3.0 x 10-7 2.5 x 10-> 4.0 x 10-° 1.2 x 10-* 


Therefor we can obtain the exact solution u(y, t) = x(y)g(t) = sin(77) exp(—t). 


Example 2. In this example, the functions ¢(y,t), ho(t), hi(t), and x(y) are 
chosen such that 


v(7,t) = (sin(ry) + cos(wy) + 2y — le~, 


with initial condition 


u(y, 0) = Aiv(7,T1) + XV) 


is the exact solution of the problem; see [22]. 

The maximum norms of the absolute error (E..) obtained with T = 
1,7, = 0.9, and 43 = —10, computed for various values of n using the 
truncated series (17) and the method in [22] are compared in Table 4. In 
addition, CPU time for different values of n is listed in this table. In Figure 
2, the exact and approximate solutions of u(7,t) in t = 0.5 for n = 6 are 
shown. 

Example 3. In the last example, we solve equations (1)—(4) with N = 1 and 
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approximation 


Out[36]= 


Figure 2: Exact and approximate solutions of u(7,t) in ¢ = 0.5 for Example 2. 


Table 4: Maximum norms of the absolute error for Example 2 


Methods in [22] Our method 
The best reported F, Fic GPU 
Time 
7 = 1/200, h = 1/400 5.946431 x 10-4 n=3 1.88002 x 10-4 0.297 s 
7 = 1/1000, h = 1/400 2.404325 x 10~° n=5 2.25173 x 10-7 ~—-0.343 s 
a 1/2000, h = 1/400 2.424209 x 107° n=7 2.27617 x 10°° 0.375 s 


1? sin(m7y) ex 
yh = oe) 


x(y) = 90, 
ho(t) = 0, 
hy(t) = 0, 


for which the exact solution is[10]: 


soe = 


The results of errors obtained with T = 1 and \,; = e~! computed for 
various values of n together with CPU time using the method presented in 
this article with the truncated series (2), with A(y) = sin(7y) are shown 
in Table 5. In Figure 3, the exact and approximate solutions of v(7,t) in 
¢ = 0.5 for n = 7 are shown. The absolute difference between the exact and 
approximate solution for n = 5 is shown in Figure 4. 
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0.6 F 4 
0.5F 4 
0.4F 4 

approximation | 
0.3 - 4 
0.2f 4 
0.1 4 
0.0 F 4 
0.0 0.2 04 0.6 0.8 1.0 
Figure 3: Exact and approximate solutions with n = 7 for Example 3. 
Table 5: Absolute errors and L? norm of errors for Example 3 
y=t n=1 n=3 n=6 n=7 
0 0 0 0 0 
0.1 4.9854 x 10-3 | 2.71955 x 107° | 1.24108 x 10-7 | 7.29989 x 10710 
0.2 9.64955 x 10-7 | 2.70053 x 10-° | 2.00455 x 10-7 | 3.9988 x 10-19 
0.3 6.75748 x 10-8 | 7.2157 x 10-® | 7.28962 x 10-9 | 4.50094 x 10-1 
0.4 1.27224 x 10-2 | 1.07862 x 10-® | 2.67932 x 10-7 | 1.05798 x 10-9 
0.5 1.29088 x 10-? | 6.40166 x 10-° | 1.47392 x 10-7 | 2.27621 x 10-1 
0.6 6.05834 x 10-3 | 5.57487 x 10-° | 2.89113 x 10-7 | 9.60188 x 10-1° 
0.7 5.56389 x 10-% | 1.56428 x 1074 | 3.7217 x 10~" | 5.04738 x 10-4 
0.8 1.61871 x 10-2 | 1.24564 x 10-7 | 1.17941 x 10-7 | 4.33895 x 10-1 
0.9 1.7427 x 10-2 | 2.72464 x 10-® | 2.59873 x 10-7 | 3.82577 x 10-22 
1 0 0 0 0 
An 6.25 x 107? 3.25521 x 10-4 | 6.78168 x 10~” | 7.56884 x 10710 
errors in L? norm | 2.11753 x 10-2 | 1.50581 x 10-4 | 3.55517 x 10-7 | 9.2949 x 10-7 
CPU Time 0.359 s 0.406s 0.437 s 0.499 s 


5 Conclusion 


In this article, the satisfier function in the Ritz-Galerkin method for a 
parabolic equation with the nonlocal initial condition has been investigated. 
Theorems on the uniqueness of product satisfier functions have been proved. 
Satisfier functions enable high adaptability to satisfy, boundary, and non- 
standard initial conditions. In addition, the use of satisfier function in the 
Ritz-Galerkin method has reduced the number of basis functions needed to 
find an accurate approximation of the problem. It is indicated that the sat- 
isfier function in the special separable case eventuated in an exact solution 
of the problem. Obtained results in numerical examples are compared with 
the results of other published methods and demonstrated the more efficiency 
of the proposed scheme. 
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_ Error 


Figure 4: The graph of the absolute error with n = 5 for Example 3. 
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